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Q ■ Abstract 

(D 

^/-v ■ We compute for the first time the decay width of A-hypernuclei in a rela- 

tivistic mean-field approximation to the Walecka model. Due to the small 
mass difference between the A-hyperon and its decay products — a nucleon 
^ \ and a pion — the mesonic component of the decay is strongly Pauli blocked 

in the nuclear medium. Thus, the in-medium decay becomes dominated by 
(3 ' the non-mesonic, or two-body, component of the decay. For this mode, the 

CS , A-hyperon decays into a nucleon and a spacelike nuclear excitation. In this 

Q"^ \ work we concentrate exclusively on the pion-like modes. By relying on the 

^|> \ analytic structure of the nucleon and pion propagators, we express the non- 

mesonic component of the decay in terms of the spin-longitudinal response 
function. This response has been constrained from precise quasielastic (p, n) 
^ , measurements done at LAMPF. We compute the spin-longitudinal response in 

r^ \ a relativistic random-phase-approximation model that reproduces accurately 

the quasielastic data. By doing so, we obtain hypernuclear decay widths 
L^ ■ that are considerably smaller — by factors of two or three — relative to existing 

iJ ■ nonrelativistic calculations. 






I. INTRODUCTION 

Understanding modifications to hadronic properties — such as masses and decay widths — 
in the nuclear medium occupies an important place in nuclear physics. The reason for 
such a status is the universal character of the problem. In electromagnetic phenomena 
one asks, for example, if the coupling of the photon to the proton gets modified in the 
nuclear environment. Similarly, hadronic processes — such as the {p,n) reaction — search for 
evidence for new states of matter through the modifications of meson (e.g., pion) properties 
in the nuclear medium. Finally, the modification of masses and decay widths for vector 
mesons might provide, through relativistic heavy-ion experiments, conclusive evidence for 
the formation of the quark-gluon plasma. 

In the present paper we are interested in understanding medium modifications to the 
properties of the lamb da-hyp eron; in particular modifications to its decay width. To our 
knowledge, this is the first time that such a study will be carried out within the framework 
of a relativistic mean-field approximation. There are several reasons why the A-decay is 
interesting. First, the mesonic mode, A —* Ntt, which comprises nearly 100% of the decay 
in free space, is unavailable in the medium because of Pauli blocking [|I|. Second, a new 
nucleon-stimulated mode, AN -^ NN, opens up — precisely because of the presence of the 
nuclear medium |0J^. As this "non-mesonic" mode dominates, the decay of the lambda in 
the nuclear medium becomes sensitive to nuclear excitations. Finally, the decay is interesting 
because the empirical A/= 1/2 rule — nearly exact for the free decay — appears to be violated 
in the nuclear medium [^]. 

Because of simplicity, or perhaps because the weak AA^vr couplings seem to be the only 
ones that are model independent |^,p, most — although certainly not all — nonrelativistic 
calculations to date have concentrated exclusively on pion-like nuclear excitations [[?|-[T0|1. In 



this, our first, relativistic study of the decay we will proceed analogously; yet our results 
seem to indicate the need for the inclusion of additional nuclear excitations p|,|6|JTT| . 

An important point that we will stress throughout this — and any future — publication is 
the need for consistency between theoretical calculations and, seemingly unrelated, exper- 
imental data. For example, pion-like modes have been studied extensively by a variety of 
experiments. Arguably the most complete set of experiments were performed at the neutron- 
time-of-flight (NTOF) facility at LAMPF []r^lT3|. These quasielastic (p, n) measurements 



placed strong constraints on the (pion-like) spin-longitudinal response, by showing conclu- 
sively that the long-sought pion-condensed state does not exist. Although improving fast, 



the present experimental situation on hypernuclear decay is still limited |T^-|T6[. Thus, 
the theoretical guidance provided by alternative experiments should not be underestimated. 
Indeed, in this work we compute the decay of A-hypernuclei using the same exact model 
that was used successfully in describing the spin-longitudinal response extracted from the 
quasielastic (p, n) experiments [IT? . 



We have organized our paper as follows. In Sec. || we develop the formalism needed to 
compute the in-medium decay of the A-hyperon. We rely strongly on the analytic structure 
of the nucleon and pion propagators to show that the non-mesonic component of the decay 
is sensitive to the spin-longitudinal response. The underlying AA^tt dynamics is prescribed 
through an effective-Lagrangian approach. For the nucleon and lambda propagators we use 
a standard mean-field approximation to the Walecka model, while the pion propagator is 



evaluated in a random-phase approximation (RPA) — using a pseudovector representation 
for the ttNN vertex. Our results are presented in Sec. |Tl|. Strong emphasis is placed on 
the dynamical quenching of the effective nNN coupling in the nuclear medium, and on the 



corresponding reduction of the hypernuclear widths. Finally, in Sec. |3 we summarize our 
findings and make suggestions for future work. 

II. FORMALISM 

The dynamics of the lambda-nucleon-pion (AA^vr) system is described in terms of the 
following effective Lagrangian Pj 

C = gwi>N{x)(l + fi;7^)r ■ 7r{x)ilJA{x) + h.c. , (1) 

where 5^^ = 2.35 x 10~^ is the weak AA^tt coupling constant and « = — 6.7 is the ratio of the 
parity conserving to the parity violating coupling. Note that — in accordance to the empiri- 
cal A/ = 1/2 rule — the lambda-hyperon field occupies the lower entry of the two-component 
"isospinor" ip\{x); the first entry is kept empty. This effective Lagrangian will now be 
used to compute modifications to the propagation of the A-hyperon in the nuclear medium. 
Specifically, the modified A-propagator in symmetric nuclear matter satisfies Dyson's equa- 
tion ig: 



Sip) = Soip) + So{pMp)Soip) + ... = Soip) + Soip)J^ip)S{p) , (2) 

where the lambda (proper) self-energy is defined by 

- ^^{P) = ^9l 1^,(1- ^l')G{p - g)(l + K^')A{q) . (3) 

Note that we have introduced the nucleon (G) and pion (A) propagators. Moreover, the 
factor of "3" in front of the integral appears as a direct consequence of the A/= 1/2 rule — the 
pTT~ channel contributes twice as much as the nii^ one. 

A. Nucleon Propagator 

The nucleon propagator will be computed in nuclear matter using a mean-field approxi- 



mation to the Walecka model |T^. Since most of the formal aspects of the derivation can be 



found in the textbook by Serot and Walecka |20] — as well as in many other publications — we 



limit ourselves to a brief review of the central issues. 

In a self-consistent (or Hartree) approximation the relativistic nucleon propagator in 
nuclear matter may be written as 



G{k) = [l/c'+M] 



N) 



where the following quantities have been introduced: 



+ T^0ik^-\k\)5{k'-E^j^\k)) 



(4) 



k'^- = (A;°-y^, k) , E^jf\k) = EU^) ± Vn , E*^{k) = ^Jk^ + Mf . (5) 

Note that once the Fermi momentuin kp is specified, the effective nucleoli mass Mj^, the 
energy shift Vat — and with them the full nucleon propagator — can be determined [^]. The 
three independent parameters of the Walecka model are the mass of the scalar meson, and 
the scalar and vector couplings respectively (the vector mass is kept fixed at its physical value 
of 783 MeV). These three parameters — which are used to describe ground-state properties 
of closed-shell nuclei throughout the periodic table — are given by: 

m, = 520MeV; gl/A'K = S.72A- gl/Aix = Ih.lhA . (6) 

The above nucleon propagator can be written alternatively in a spectral representation. 
That is, G{k) = G^+\k) + G^-\k), where 

'^'"^''^ V*(k,s)ir(k,s), (7a) 




-ir] k^ - E)^> {k) + ir] 
k, s 



'jv; -r t/// 

(7b) 

k^ + E'}^>{k)-i'q 

We have chosen to display the analytic structure of the nucleon propagator to highlight the 
relevant physics that has been included at the mean-field level. For example, the G^~^ term 
represents the anti-particle (or filled Dirac sea) contribution to the nucleon propagator. This 
term is analytic in the full complex fc*^-plane, except for the presence of negative-energy poles 
located (slightly) above the real axis. This component of the propagator does not contribute 
to the width of the A-hyperon in the nuclear medium. The G'^^^ term, in contrast, represents 
the positive-energy part of the propagator. In addition to a change in the dispersion relation 
because of the presence of the mean fields, the only other modification to this part of the 
propagator, relative to its free-space structure, is the shift in the position of the pole of all 
the states below the Fermi level from slightly below to slightly above the real A;°-axis (see 
Fig.|l]). In this way, the "quasinucleon" propagator takes a form analogous to the one in 
free space — the conventional Feynman propagator — while at the same time enforcing the 
Pauli-exclusion principle. 

B. Pion Propagator I. 

Isolating the analytic structure of the pion propagator in nuclear matter is not as straight- 
forward as in the case of the above nucleon propagator. Thus, in this case it is convenient 
to express the pion propagator in its more general form — by using its Lehmann representa- 
tion ig: 

7r JO oj — q^ — iTj IT J-oo oj — q^ + iTj 

In this form, knowledge of the spectral weight (i.e., the imaginary part of the propagator) is 
sufficient to reconstruct the full propagator, albeit at the expense of introducing an additional 
integral. Yet the virtue of such a representation is that the physical content of the imaginary 
part is simple and illuminating. We will return to this point later in the section. Now we 
proceed to compute directly the imaginary part of the lambda self-energy. 



C. Lambda Self-Energy 

As the lambda propagates through the nuclear medium it feels the presence of the strong 
scalar and vector mean fields. Thus, the lambda propagator gets modified at the mean-field 
level in a manner similar to that of the nucleon propagator. However, because there are no 
pre-existing lambdas in the nucleus, there is no Pauli correction to the propagator. Hence, 
the lambda propagator looks exactly as a conventional Feynman propagator 

but with a mass and a dispersion relation modified by the strong mean fields: 

pM^^(pO_^^^p)^ E^^\p)^EXip)±VA, EXip)^^p^ + Mf. (10) 

The only two parameters that remain to be specified are the coupling constants of the 
lambda to the scalar and vector fields. Here we determine them from the assumption that 
the strange quark in the lambda does not couple to the sigma nor to the omega meson. In 
this way we end up with the simple quark-model estimates of: 

QcrKA/gcrNN = guj\\/gujNN = 2/3 . (11) 

Of course other possible values — perhaps phenomenologically more robust — can be easily 
incorporated in our calculation. Yet, irrespective of these values the effective lambda pa- 
rameters become: 

(MI-M^) = (^^{M;,-M^); Va=(^^V,,. (12) 

\9crNN / \gujNN / 

In a mean-field approximation it is the strong scalar field that is responsible for inducing 
a shift in the mass of the lambda, relative to its free-space value. Yet, its entire decay 
width is generated from the weak matrix element given in Eq.(^; note that the dispersive 
contribution to the mass term is insignificant. Moreover, the contribution to the width 
comes entirely from the particle i+if]) component of the nucleon propagator [the last term 
in Eq.(|7a])]. This in turn constrains the contribution from the pion propagator to arise 
exclusively from the first {—irj) term in Eq.(|]); otherwise the g° integral in Eq.(^ would 
vanish. In this way, the q^ integral can be computed with the aid of Cauchy's theorem. We 
obtain, 

X.S(p) - 3g^J ^^ ^^^ e[E^{k) - E_ 

J^^dujI,^A{q,uj)6[E^^\p)~Ei;\k)-u) , (13) 

where now k^* is evaluated at its on-shell value: k^* = {E^(k),'k = p — q). To proceed, and 
in particular to isolate the Lorentz structure of the self-energy, it is convenient to use the 
following identity between gamma matrices: 

{1-K^')i^* + M*^){l + K^') = (1-«:2)M^ + {1 + K^)lt'' + 2fi:|^V • (14) 



Note that a parity-violating — axial-vector — component has been generated at the one-loop 
level. Although interesting, the effect of this term to the decay will be of 0{g^) and will be 
neglected henceforth. In this way we arrive at the following form of the A-propagator: 



Sip) 



Ao7° - Av7 ■ P + 5 



^0 " ^v 



fi2 



where the following quantities have been introduced 



.0* 



^v — IpI ^ ^v ; 

B = M 



A 



(15) 



(16a) 
(16b) 
(16c) 



All these quantities depend (for an on-shell lambda) on the momentum of the lambda (p) 
and the density of the system (kp); for simplicity, reference to these parameters has been 
suppressed. The medium-modified width of the A-hyperon is now extracted by evaluating 
the imaginary part of the denominator in Eq.(^) at the position of the pole. That is. 
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(17) 



Since in our model the medium-modified pion propagator depends solely on the magnitude 
and not on the direction of q, the delta function in Eq.(|l^) serves to perform the angular 
integrations. This results in the following form for the imaginary part of the A-self-energy: 



TnXiiji) 
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where the Cj-coefficients are defined by 



Ci{q,uj) 



'1- 
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-K 



^^.(p-q) , 

p| -qp 



--I^A{q,uj) 
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for z = S ; 
for z = ; 
for i = V , 



(18) 



(19) 



supplemented with the energy-momentum-conservation relation 

E]V(p-q) = /(p^^qm¥ = ^^^^(p-q)-^^=4^^(p)-^-v]v. 



(20) 



Note that suitable kinematical conditions at the various interaction vertices impose con- 
strains on the limits of integration for both the q- and cj- integrals in Eq. ([T8|) . For a detailed 
discussion on these limits we refer the reader to the appendix. 

D. Pion Propagator II. 

Equation (|18|) is very general, as it depends solely on the analytic structure of the pion 
propagator. In this section we discuss in detail the relativistic model used to calculate the 
pion propagator — placing special emphasis on the physics underlying its imaginary part. 
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The imaginary part of the pion self-energy is interesting and fundamental as it is related, 
after a suitable analysis, to a physical process: the (plane-wave) cross section in proton- 
neutron scattering. There are a variety of physical processes that modify the propagation 
of the pion as it moves through the many-body environment. These include the coupling of 
the pion to particle-hole (ph) and nucleon-antinucleon {NN) excitations. In nuclear matter 
these two kind of excitation have a distinctive character. For example, pions with spacelike 
(g^<0) momenta can only "decay" into (real) ph-pairs. In contrast, the A^iV-channel opens 
up at the relatively large timelike momentum of q'^ = 4M^>0. 

As was discussed in the Introduction, the mesonic decay of the lambda — in which an 
on-shell pion with timelike momentum is produced — is strongly Pauli-suppressed in nuclear 
matter [this process is represented by diagram (b) in Fig. |]. Recall that in free space the 
nucleon is produced with a momentum of only 100 MeV, which is well below the Fermi 
momentum at saturation. However, with the closing of the mesonic branch the non-mesonic 
channel opens up. For this nucleon- induced (AN^NN) decay the exchanged pion is con- 
strained to have spacelike momentum. Thus, the coupling of the pion to spacelike particle- 
hole excitations becomes of paramount importance [this channel is represented by diagram 
(c) in Fig. 1^. We now compute the medium- modified pion propagator in nuclear matter 
by iterating — to infinite order — the lowest-order contribution to the self-energy by means of 



Dyson's equation ||T8i. That is. 



A(g) = Ao(g) + Ao(g)no(g)A(g) = Ao(g) + Ao(g)n(g)Ao(g) , (21a) 

n(g) = no(g) + no(g) Ao(g)n(g) = no(g) + no(g)A(g)no(g) . (21b) 

where we have introduced the free pion propagator 

Ao(g) = ^j— - , (22) 



and the lowest-order pion self-energy 



The factor of A = 2 represents the isospin degeneracy in symmetric nuclear matter and 
f^/4:7i = 0.0778 represents the strength of the pseudovector nNN vertex. Further, we have 
established the equivalence between various ways of computing the medium-modified pion 
propagator. It is the last form in Eq. (|21a|) that we use here to write the mesonic and 



non-mesonic contribution to the lambda self-energy: 

Im^iip) = -^-ri / QdQ dujCi{q,uj) 



--ImAo{q,uj] 



(24a) 



J^S,(p)^^^^ =-^r^\dql duQiq,uj)Aliq,uj)SLiq,uj). (24b) 



Note that we have introduced the spin-longitudinal response Si{q,u!), to be defined below. 
The imaginary part of the lambda self-energy is obtained by "cutting" all the Feynman 
diagrams given in Fig. |l], or equivalently, by putting all intermediate particles on their mass 



shell. All of these diagrams contain a nucleon in the intermediate state that must lie above 
the Fermi surface. Certain diagrams [such as diagram (b)] have a pion in the intermediate 
state and contribute to the mesonic component of the decay. However, as the mesonic 
component is known to be strongly suppressed by the Pauli exclusion principle, no effort 
has been made to compute modification to the pion mass in the nuclear medium [note that 
even in the chiral {m.,^ = 0) limit the momentum of the outgoing nucleon will increase to only 
160 MeV; still well below the Fermi momentum at saturation]. The decay of the lambda into 
a nucleon and a (spacelike) particle-hole excitation [such as in diagram (c)] constitutes the 
non-mesonic component of the decay. In this manner, the non-mesonic component becomes 
proportional to the imaginary part of the lowest-order pion self-energy, or equivalently, to 
the spin-longitudinal response: 

S^^\q,uj) = --Z^Uo{q,uj). (25) 

Finally, diagrams such as in Fig. |l](d) contribute to the non-mesonic decay through the 
higher-order iteration of the lowest-order pion self-energy. These diagrams have been 
summed to all orders using an RPA approach. Note that for the residual interaction we 



have included a phenomenological Landau-Migdal term, as described in detail in Ref. [ff7 



III. RESULTS 

This section is devoted to the discussion of our results. Before doing so, however, we 
clarify the various approximations required to extend our nuclear-matter results to finite 
nuclei. 

The first step in such a procedure is to calculate the ground state of various closed-shell 
nuclei (such as ^''Ca) in a relativistic mean- field approximation to the Walecka model p!9| , p0 



Once the self-consistent mean fields are generated, we solve the Dirac equation for a A- 
hyperon moving in these potentials; no additional self-consistency is demanded. Hence, the 
lambda moves in a mean-field potential that differs from the corresponding nucleon potential 
in only two ways. First, the scalar and vector potentials were scaled by a factor of 2/3, as 
the strange quark in the lambda is assumed not to couple to the sigma nor to the omega 
meson. Second, by assuming that the magnetic moment of the lambda is carried exclusively 
by the strange quark, a strong tensor mean-field potential is generated. To illustrate these 
modifications, we have included in Table | a comparison between the nucleon and the lambda 
single-particle spectra using a mean-field potential appropriate for ^°Ca. For the case of the 
s^^^ states — which lack spin-orbit partners — the binding energy for the lambda orbitals is 
almost 2/3 that of the nucleon. Yet the most interesting modification is the near total 
dilution of spin-orbit effects in the lambda spectrum due to the strong tensor interaction. 

Once the A-spectrum has been obtained, we compute the effective (or average) vector 
density sampled by a A-hyperon occupying the lowest (s^/^) orbital. That is, 

(pv) = J rf V(r)f7A(r)70f/A(r) . (26) 

This is then the nuclear-matter density that is used to compute the decay width: rA(p; (Pv))- 
Finally, the A-hypernucleus decay width is extracted from the convolution of this expression 
with the momentum distribution of the lowest orbital: 
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/^rA(p; (p.))f/A(p)rf/A(p) . (27) 



A. Mesonic Decay 

Computing the mesonic component of the decay is relatively straightforward. We start 
by writing the free pion propagator of Eq.(p2[) in its spectral form 



Ao(g) 



2t^q 



uj^^Ja^ + ml. (28) 



00 — UJq + trj 00 + UJq — ZT] 

This yields the following simple form for the mesonic component of the self-energy [Eq.( p4a| )] 

Im^Ap) = — ^^ r dooCAoo) . (29) 

Here the constant Cj introduced earlier must be evaluated at its "on-shell" value oo = oOg. 
Moreover, the limits of integration are obtained from simple kinematical considerations after 
using the definitions given in Eqs. (|A2|) and ( |A3| ) . As the weak AiVvr coupling constants have 



been chosen to reproduce the width of the lambda-hyperon in free space, we obtain in the 
limit of /cf = 0: 

To = 2.490 X 10"^^ MeV . (30) 

Yet the mesonic component of the decay gets modified dramatically in nuclear matter be- 
cause of strong Pauli correlations. This can be seen in Table |l| which shows that the mesonic 
component of the decay drops dramatically with baryon number. We will return to discuss 
these results in more detail after presenting results for the non-mesonic component of the 
decay. 

B. Non-Mesonic Decay 

The non-mesonic component of the decay is considerably more complicated to evaluate. 
Not only does one need to compute the ground-state of uniform nuclear matter but, in ad- 
dition, one must compute its linear (spin-longitudinal) response. The linear response of the 
ground state — contained in the imaginary part of the polarization insertion — requires a spe- 
cific representation for the irNN vertex. In this work we have adopted a pseudovector (PV) 
representation. The pseudovector representation adopted here is just one of several possible 
choices. Indeed, we could have adopted instead a pseudoscalar (PS) representation, which is 
guaranteed to be equivalent on-shell. The merit of the PV representation, however, is that 
the correct soft-pion limit of various observables is enforced at the level of the Lagrangian 



density, rather than from a sensitive cancellation among various Feynman diagrams [ pi[] . 

One of the most interesting results that emerged from our earlier work on the pion- 
nucleon system is the dynamical quenching of the nNN coupling constant in the nu- 



clear medium |T^J2^. The origin of this quenching is intimately related to the the non- 
conservation of the nucleon axial- vector current or, equivalently, to the presence of a nucleon- 
mass term. That is, 
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q,Uip')ri'Uip) = 2MMU{p'h'U{p) ^ . (31) 

In the nuclear medium the nucleon mass is reduced — dynamically — from its free-space value 
because of the strong scalar field. Thus, relative to a "nonrelativistic" calculation, by which 
we mean that the potentials are small and M^ -^ Mn, it appears as if the ttNN coupling 
constant has been effectively reduced in the medium. Moreover, this reduction is density 
dependent and tracks the behavior of the effective nucleon mass; namely, the quenching 
increases with the density of the system. This relativistic behavior was responsible for gen- 
erating no pion condensation — even in the absence of short-range correlations [^. More 
importantly, it described correctly [|l^ the behavior of various quasielastic {p, n) spin observ- 



ables measured at LAMPF [p!2| , p!3[ . Measuring a variety of spin observables was required for 
the extraction of the the spin-longitudinal response, which was instrumental in eliminating 
pion condensation as a possible new state of matter. With these results in mind, we now 
proceed to use the same exact relativistic model to compute the non-mesonic decay of the 
A-hyperon. 



Equation (|24b| ) suggests that once the spin-longitudinal response is obtained, computing 



the different Lorentz components of the non-mesonic decay becomes straightforward. We 
present in Fig. ^ results for the decay width of the A-hyperon at a nuclear-matter density 
{kp = kp) appropriate for ^^Ca (see discussion at the beginning of the section). The momen- 
tum distribution for the lowest (s^^^) A-orbital is displayed by the solid line in the figure; 
it is normalized so that the area under the curve — which has been multiplied by a factor of 
100 for clarity — is equal to one. The other three curves in the figure display the mesonic and 
non-mesonic components of the decay computed in an approach that has fixed all masses to 
their free-space value; this is our best attempt at reproducing nonrelativistic results. The 
dashed curve shows the mesonic component of the decay as a function of the momentum of 
the lambda. Clearly, the mesonic mode is strongly suppressed for momenta below 200 MeV. 
The curve eventually "heals" to its free-space value, but too late to make a contribution 
to the in-medium decay; we have computed a mesonic contribution to the decay in ^°Ca, 
relative to its free-space value, of only 0.014. This number was obtained by folding the de- 
cay width computed in {M^ = M]\f) nuclear matter with the momentum distribution of the 
s^' ^-orbital (see inset). The non-mesonic contributions to the decay, computed without and 
with RPA correlations, are displayed by the dot-dashed and dotted lines, respectively. These 
contributions peak at small momenta and, thus, have good overlaps with the momentum 
distribution. The Hartree (or non-RPA) and RPA contributions to the width have increased 
to 0.834 and 1.053, respectively. 

The same calculation has been repeated in Fig. ^, but now with self-consistent nucleon 
and lambda propagators. That is, in-medium masses were used with their values determined 



self-consistently from the mean-field equations |20]. The different scaling of the nucleon and 
lambda masses with density generates a slight change in the mesonic component of the decay, 
but one that it still too small to impact the total decay width. By far the most important 
modification emerges from the dynamical quenching of the nNN coupling in the nuclear 
medium [^ . Indeed, the non-mesonic components of the decay computed in the relativistic 
approach are now reduced to 0.444 and 0.544, respectively. This large reduction factor, 
relative to the above M^ = M^ calculation, is consistent with the reduction observed in the 



quasielastic (p, n) reaction [0. Moreover, the relativistic RPA result obtained here for the 



total width — F/Fq = 0.567 — is more than a factor of three smaller than the nonrelativistic 
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value reported by Ramos, Oset, and Salcedo in Ref. 0. This large reduction factor is not 
exclusive to ^°Ca, but is seen in all hypernuclei reported in Table |l[ 

It is also relevant to mention that in a very recent publication (of only a few days old!) a 
new nonrelativistic evaluation of the A-decay width has been made [^. For the case of ^'^Ca 
a total decay width (without including the two-particle — two-hole component of the decay) 
of r/Fo = 1.08 has been reported. This represents a reduction of about 35-40% relative 
to the value presented in Ref. [0. Yet, it is essentially identical to our "nonrelativistic" 
(M^ = Mtv) value of F/Fq = 1.07. However, had we used the same M^ = Mj^ model to study 



the quasielastic (p, n) observables IT^JTH] , we would have grossly overestimated the observed 



spin-longitudinal response [T^ 



IV. CONCLUSIONS 

We have computed for the first time the decay width of A-hypernuclei in a relativistic 
mean-field approximation to the Walecka model. Most of the formalism presented here 
is quite general, as it relies only on the analytic structure of various propagators. The 
calculation was performed in two stages. First, ground-state properties of various closed- 
shell nuclei were computed self-consistently. With the mean-fields potentials in hand, the 
Dirac equation for a A-hyperon was solved with scaled scalar and vector potentials. A scaling 
factor of 2/3 was introduced from the assumption that the strange quark in the lambda does 
not to couple to the sigma nor to the omega meson. Moreover, by assuming that the magnetic 
moment of the A-hyperon is carried solely by the strange quark, a tensor mean-field potential 
was introduced that diluted — almost completely — all spin-orbit effects. Second, the mesonic 
and non-mesonic components of the decay width — with and without RPA correlations — were 
computed in infinite nuclear matter. The uniform nuclear-matter density was estimated from 
evaluating the average density sampled by a A-hyperon occupying the lowest (s^/^) orbital. 
Finally, the nuclear-matter results were folded with the momentum distribution of the lowest 
A-orbital to compute the finite-hypernuclear width. 

As expected, the mesonic component of the decay was strongly suppressed by Pauli cor- 
relations. Indeed, for the values of the nuclear-matter density adopted here, the decay of 
A-hyperons with momenta less than 200 MeV was strictly forbidden. Since in the Walecka 
mean-field model the momentum components of the A-hyperon beyond 200 MeV are small, 
the strong suppression of the mesonic mode ensued. In contrast, the non-mesonic component 
of the decay is not Pauli-blocked. In the non-mesonic mode the A-hyperon decays into a 
nucleon and a pion-like pair with spacelike momenta. Information on these kind of pion-like 
excitations is contained in the nuclear spin-longitudinal response. This fundamental nu- 
clear response is accessible through charge-exchange {p, n) reactions. Indeed, high-precision 
measurements made at the NTOF facility at LAMPF were instrumental in answering some 
fundamental nuclear-physics questions |T2|,|T3|. For example, these experiments showed con- 
clusively, by extracting the spin-longitudinal response, that the long-sought pion-condensed 
state does not exist. In our relativistic model the absence of pion condensation was attributed 
to a dynamical quenching of the nNN coupling constant in the nuclear medium |]17|^ . 

In the present work we have computed the non-mesonic component of the decay — using 
the same relativistic RPA model that was used to confirm the LAMPF measurements. The 
dynamical quenching of the nNN coupling constant in the medium was now responsible 
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for generating decay widths for A-hypernuclei that were two to three times smaller than 
those obtained by existing nonrelativistic calculations |9.23] — and considerably smaller than 



measured experimentally |T^. This suggests that other modes of excitation, besides pion-like 



modes, might be important in understanding the decay width of A-hypernuclei P,§,^ 



The study of such modes seems both interesting and necessary. For example, a weak 
AA'^o; vertex will induce a coupling to isoscalar-vector modes. As in the pion case, this 
additional contribution to the non-mesonic width can be related to a well-known nuclear 
response: the longitudinal response measured in electron scattering. Hence, its impact on 
the decay will be also strongly constrained. Searching for the underlying dynamics behind 
the decay width of A-hypernuclei — while at the same time maintaining consistency with the 
wealth of existing nuclear-response data — will be the subject of a future report. 

ACKNOWLEDGMENTS 

This work was supported by the DOE under Contracts Nos. DE-FC05-85ER250000, 
DE-FG05-92ER40750. 



12 



APPENDIX A: KINEMATICAL CONSTRAINTS AT THE VERTICES 

The limits of integration in the q- and cj- integrals of Eq. (|^) are constrained by energy- 
momentum conservation at the various interacting vertices. We now examine the kinematical 
constraints imposed at each vertex separately. 

1. Kinematical Constraints at the Lambda-Nucleon-Pion Vertex 

Energy-momentum conservation at the AA^vr vertex demands that an on-shell lambda- 
hyperon "decays" into an on-shell nucleon — above the Fermi surface — and, in principle, an 
off-shell pion. That is, 

^ = 4+Hp)-E^+Hp-q)>o. (Ai) 

Moreover, as the angular integration over q can always be performed, the following limits 
of integration ensue: 



'^max 



4+)-max(E^+)(A;p),Ei+)(p + g)) , (A2a) 

4+)-max(E^+)(fcp),E^+)(p-g)) , (A2b) 



where E)^ '{p ± q) = J{\p\ ± q^ + M^ + Vn- Note that the above relations impose con- 
straints on the g-integral as well. For example, the w-integral will vanish unless ujma.x> ^min- 
Further, u^i^ is constrained to be greater than or equal to zero. Thus, from these two 
conditions the following limits are obtained: 



n(+) T/..^2 /\/r*2 



gi„in = max(0,A;F-|p|) , gmax = max 0, V(^a -'^Jv)^ - M]^2 _ |p| \ ^^3) 



2. Kinematical Constraints at the Particle-Hole-Pion Vertex 

In the case of the non-mesonic decay of the lamb da-hyp eron, additional constraints — at 
least in nuclear matter — follow from enforcing energy-momentum conservation at the phyr 
vertex. Indeed, the energy of the exchanged (spacelike) pion is given by: 

u = E^+) (k + q) - E^+) (k) > , ( A4) 

where E}^ (k + q) > E}^ (kp) denotes the energy of the particle and E}^ (k) < E}^ {kp) is 
the energy of the hole. As in the above case, the most stringent constraints emerge in the 
case in which the particle, the hole, and the pion are coUinear. We obtain. 



tfn 



0J„ 



max (^0, ^(|q| -A;^)2 + M*^ - ^kj, + M*^^ , (A5a) 



{H + kpf + M^^ - Jkl + Mf . (A5b) 



The resulting limits of integration for the cu- integral, then follow from enforcing Eqs. ([A^ 
and (|A5|) simultaneously. 
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FIG. 1. Spectral content of the nucleon propagator in a relativistic mean-field approximation. 
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FIG. 2. Feynman diagrams contributing to the in-medium decay width of the A-hyperon. 
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FIG. 3. Momentum dependence of the A-hyperon decay width, normalized to its free-space 
value. The mesonic contribution to the width is displayed with the dashed line, while the 
non-mesonic decay — with and without RPA correlations — is contained in the dotted and dot-dashed 
curves, respectively. All calculations were done assuming the "nonrelativistic" {M^ ^ Mi\f) limit. 
The inset shows the same curves but now folded with the momentum distribution of the lowest 
(s^'^) lambda orbital in ^^Ca. 
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FIG. 4. Momentum dependence of the A-hyperon decay width, normalized to its free-space 
value. The mesonic contribution to the width is displayed with the dashed line, while the 
non-mesonic decay — with and without RPA correlations — is contained in the dotted and dot-dashed 
curves, respectively. All calculations were done assuming the self-consistent (M^) value for the 
effective nucleon mass. The inset shows the same curves but now folded with the momentum 
distribution of the lowest (s^'^) lambda orbital in ^^Ca. 
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TABLES 

TABLE L Single-particle binding energies (in MeV) for a nucleon or a lambda in the 
self-consistent mean-field of ^'^Ca. Quark-model estimates were used for go-AA and QojAA- See 
text for details. 



Orbital 


N 


A 


,1/2 


-55.390 


-36.274 


p3/2 


-38.903 


-24.030 


pl/2 


-33.181 


-23.591 


d'/'' 


-22.751 


-11.769 


,1/2 


-14.389 


- 9.838 


d3/2 


-13.876 


-11.099 



TABLE n. Hypernuclear decay widths in a relativistic mean-field approximation to the Walecka 
model. The third column contains the mesonic contribution to the decay width, while the fourth 
and fifth columns display the dominant non-mesonic and total decay widths, respectively. The 
last two columns contain the nonrelativistic results, without including the two-particle — two- hole 
component of the decay, as computed in Ref. M and in Ref. p^, respectively. 



Nucleus (fcF)A^ Tm/Tq T^M./Tq TTotal/rp T/Tq g T/Tq g 

1.07 
1.09 
1.08 



^^c 


0.991 


0.112 


0.413 


0.525 


1.76 


\^0 


0.956 


0.099 


0.497 


0.596 


1.78 


2«Si 


1.013 


0.053 


0.479 


0.532 


— 


^A^S 


1.024 


0.048 


0.491 


0.539 


— 


\°Ca 


1.000 


0.023 


0.544 


0.567 


1.79 



